Order and Creep in Flux Lattices and CDWs Pinned by Planar Defects 
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The influence of randomly distributed point impurities and planar defects on the order and 
transport in type-II superconductors and related systems is considered theoretically. For random 
planar defects of identical orientation the flux line lattice exhibits a new glassy phase with diverging 
shear and tilt modulus, a transverse Meissner effect, large sample to sample fluctuations of the 
susceptibility and an exponential decay of translational long range order. The flux creep resistivity 
for currents J parallel to the defects is p{J) ~ exp— (Jo/J)'' with p = 3/2. Strong disorder enforces 
an array of dislocations to relax shear strain. 

PACS numbers: 74.25.Qt, 71.55.Jv, 74.62. Dh, 64.70.Rh 



Introduction. Type-II superconductors can be pene- 
trated by an external magnetic field in the form of quan- 
tized magnetic flux lines (FL). Under the influence of a 
transport current J FLs will move and hence give rise 
to dissipation. The resulting linear resistivity is propor- 
tional to the magnetic induction B To stabilize su- 
perconductivity it is therefore essential to pin FLs. One 
source of pinning is point disorder. In high-Tc materi- 
als point disorder is practically always existing because 
of the non-stoichiometric composition of most materials. 
Then the system regains superconductivity in the sense 
that the linear resistivity vanishes 2]. However, ther- 
mal fluctuations lead to flux creep resulting in a non-zero 
nonlinear resistivity of the form p{J) ~ e"^'^^/'^-'^ where 
^ = 1/2 Q . Jp{'> J) is a function of i3, temperature T 
and the concentration and strength of the pinning cen- 
ters. This response of the system on an external current 
is closely related to the order of the FLL in the presence 
of point pinning centers, which shows a power law decay 
of its translational order parameter in the " Bragg glass" 
phase [3, 0, 0] . Although the linear conductivity is now 
zero, there is still a finite resistivity for a finite current. It 
is therefore indicated to look for a more effective pinning 
mechanisms corresponding to larger values of the creep 
exponent /z. One option is columnar defects which lead 
to a " Bose glass" phase with stronger pinning properties 
(see e.g. Q). 




FIG. 1: Planar crystalographic defects in BSCCO (bright ver- 
tical lines), unpublished. Figure by courtesy of M. Menghini, 
Y. Fasano and F. de la Cruz from Centro Atomico Bariloche, 
Argentina and Eli Zeldov from Weizmann Institute, Israel. 



An even more pronounced effect can be expected from 
planar defects like twin planes or grain boundaries, which 
will be considered in the present paper. Twins are ubiqui- 
tous in superconducting YBCO and La2Cu04 where they 
are needed to accomodate strains arising from tetragonal 
to rhombic transformations. But also other causes are 
possible (see Fig.l). Planar defects occur frequently in 
families with the same orientation but random distances 
0, 0] or in orthogonal families of lamella (" colonies" ) . 
The mean distance Id of the defect planes is of the order 
of 10 nm to [lo| . Pinning of individual FLs by 
planar defects has been investigated in the past both for 
clean and disordered systems 0, [ll|, [13] ■ Recently it was 
shown that, depending on the mutual orientation of the 
FL lattice (FLL) and the defects, dilute planar defects 
are indeed a relevant perturbation even in the presence 
of point disorder [isj, provided they are parallel to the 
main lattice planes of the FLL. In systems with parallel 
defect planes this is the generic situation since the FLL 
will rotate in such a position to reach maximum overlap 
with the defects (provided B is aligned with the defect 
planes). It turns out that the Bragg and the Bose 
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FIG. 2: Disordered vortex lattices resulting from impurities, 
columnar and planar defects of concentration nimp,ncd and 
Upd, respectively. In the presence of planar defects the planar 
glass phase is ultimately stable. ^ denotes the creep exponent. 

glass phase are unstable with respect to the presence of 
many random planar defects and will be substituted by 
a new type of a planar glass phase (see Fig. 2). In the 
present paper we will discuss the nature of this phase. 
The latter is characterized by a complex energy land- 
scape with many metastable states and diverging energy 
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barriers leading to a new creep law with fi = 3/2, large 
sample to sample fluctuations of the magnetic suscepti- 
bility, an exponential supression of translational order in 
the direction perpendicular to the defects, a resistance 
against shear deformations as well as the occurrence of 
a transverse Meissner effect. If only displacements per- 
pendicular to the defects are considered, as in the main 
part of this paper, our results apply also to a wide class 
of systems which exhibit regular lattices of domain walls 
like magnetic slabs, charge density waves 13 and incom- 
mensurate systems 
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Model. We consider an Abrikosov FLL in the presence 
of randomly distributed point impurities and random de- 
fect planes, aligned with the magnetic field. Since in both 
types of imperfections superconductivity is suppressed 
they will attract FLs. Then the Hamiltonian reads [l^ 



n = 



a 

+2[yp(r) + Vb(r)]p(u,r)} 



(1) 



where a, /?, 7, 5 run over x, y. u(r) = {ux, Uy) denotes the 
FL displacement. Only components of the elastic con- 
stants Cafs-yS with pairwise equal indices are non-zero 17 1 . 



-u)} 



Pp is the 



/o(u,r) = po {-V±u 
FLL density with po — Bo/(t>Oi (f'o is the flux quantum. G 
is a reciprocal lattice vector of the FLL and r_L = ix,y). 
Vp(r) denotes the pinning potential resulting from ran- 
domly distributed point impurities. We will first con- 
sider the (realistic) case that all defect planes have the 
same orientation but random distances Q • Then the FLL 
will orient itself such that its main lattice planes will be 
parallel to the planar defects to allow for their maximal 
overlap [l3|. The defect pinning potentials then has the 
form Voi'r) = —voj^d^i^ ~ ^d) where we assumed 
that the defect planes are parallel to the yz— plane. The 
(5-functions are considered to have a finite width of the 
order of the superconductor coherence length ^c- A rough 
estimate for the defect strength is given by vd ~ ^^cCc' 
He is the thermodynamic critical field. The statistical 
properties of the pinning energies are then encoded in 
their pair correlation functions Rp{u) and Rd{ux) for 
point disorder and planar defects, respectively. Since 
the FLL density includes a slowly varying and a peri- 
odic part, ps and pp, respectively, we decompose the 
pinning energy density accordingly. From the periodic 

part we get Rd{ux) = (wdPo)^/^d E«#o e"^''''"^^ " is 
integer [l8|. £ £d, where £ and £d are the mean spac- 
ing between the FLs and the defect planes, respectively. 
The contributions from ps do not contribute to the glassy 
properties of the system since they can be eliminated by 
a simple transformation [l9| . 

Since our main concern are the defect planes, it seems 
to be justified to start with a simplified model in which 
only the displacements Ux = u of the FLs perpendicular 
to the defect planes are considered. Then only the elastic 
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C44 remain in the Hamiltonian. From a technical 
point of view it is convenient to consider a generalization 
of our model in d dimensions by replacing x hy a {d— 2)- 
dimcnsional vector x. 

Weak disorder . In the absence of defect planes point 
impurities are relevant in less than 4 dimensions. The 
FLL exhibits a phase with quasi long range order : the 
Bragg glass 3, 4,^], which exhibits a power law decay of 
5-^(1.) = (e«G(u(r)-u(o))^ _ Irj-^'^-'^). The Fourier trans- 
form of 5*0 (r) is the structure factor which has Bragg 
peaks. 

It was recently shown in [l3| that dilute planar defects 
can be a relevant perturbation also in the presence of 
point disorder. Indeed, distorting the initially ordered 
FLL in volume L'^^^LyL^ the energy gain is of the order 
— {R'l^' {0)L'^^^y^^ LzLy whereas the elastic energy loss is 
ciiLzLyL'^^* since distortions are aligned parallel to the 
defects. For L:^ Ld (c?i/i?^"(0))i/(6-'') the pinning 
energy gain wins and the FLL starts to disorder in the 
directions perpendicular to the defects. The critical di- 
mension below which weak planar defects are relevant is 
d = 6. 

For a more detailed study we use now a functional 
renormalization group approach in d = 6 — e dimensions. 
We follow closely a related approach for columnar dis- 
order [2^, [2lj but keep the unrescaled quantities which 
correspond to the effective parameters measured on scale 
L. To lowest order the flow equations for e ^ 1 read 

dlncij/rflnL = 2i^^"(0)L7(47^Cll)^ i = 4, 6 (2) 
dRD{u)/d\nL = R'Ij{u)L'' {R'o{u) - 2i?^(0))/(47rcii)^ 

Thermal fluctuations and point disorder are irrelevant 
for e < 4 and e < 2, respectively. There is no renor- 
malization of Cll because of a statistical tilt symmetry 
[2^ . For L — > L/), many metastable states appear and 
R'[, (0) develops a slope discontinuity at the origin which 
results in diverging elastic constants C44 and cee- The 
renormalization can however be continued to L ^ L]j if 
one imposes a small but finite tilt of the FLL such that 
i?""(0) has to be replaced by R""{0+) in Eq. In 
this case C44 and cqq remain finite but new terms of the 
form f d(l) l^y cos (p (dyu) + T,z sin (j){dzu)\ /M are gen- 
erated in the energy density which dominate the energy 
for small u. The fixed point function R*jj {u,L)L'^ — 

(27rcii)^e 1^ — ^ (u — 1)^ for < m < ^ is periodic in 
u with period £. The newly generated terms rcnormalize 
according to 



CQQ^'^dJ^y/dliiL ; 



c^^^'^dJ^z/dlnL 



eVcn^VlSi- (3) 



^z(y) has the meaning of a interface tension of a domain 
wall parallel to the x and y (z) axes. Y.^ can be measured 
by changing the external magnetic field by HxX which 
changes the Hamiltonian by —{Bq/Att) J d^r HxdzU. To 
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tilt the flux lines with respect to the z-axis, has to 
overcome the interface energy ~ which results in a 
threshold field H^ ,, = ^t^'^ I [i'oV^) below which FLs 
remain locked parallel to the planes. This is the trans- 
verse Meissner effect: a weak transverse magnetic field 
Hx is screened from the sample. In this case C44 is in- 
finite! Only for > H^^c the average tilt of the FLs 
becomes non-zero and C44 stays finite. Moreover, there 
is a resistance against shear of the FLL: the shear defor- 



moved by a distance £. This droplet is a saddle point, as 
usual in nucleation phenomena. In the presence of planar 
defects the energy of the nucleus has the from 
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(4) 



Here we have taken into account that the elastic energy 
and the energy from the disorder scale in the same way. 



The saddle point Ly/Y,y — Lzj^x ^ i /cn 



mation dyU^ is non zero (and cee finite) only if the shear for the nonlinear resistivity in d = 3 (J < Jd) 



/ ^ gives 



stress Oxy is larger then a critical value "^^y/i, otherwise 
C66 is infinite. The divergence of cee is a new property 
which does not exist in Bose glass. 

An infinitesimal change 5Hzi in the longitudinal field 
allows to measure the longitudinal susceptibility x — 
B[)d{dxu) / dHz- The disorder averaged susceptibility 
X = Bq/Attch is independent of the disorder as a result of 
the statistical tilt symmetry. The glassy properties of the 
systems can most easily be seen by the sample to sample 
fluctuations of the magnetic susceptibility x^ ~ X^- Per- 
turbation theory gives (x^-X^)/x^ = i?^"(0)L7(5c?i) - 
{L/LuY, i.e. the sample to sample fluctuations of the 
susceptibility grow with the the scale L < Lu, d < 6 
which is a signature of a glassy phase [l^l . 

The structural correlations in this phase are obtained 
in the standard way from R'ly{u,L) [| which gives 
S'g(x, y, z) ^ Ixl^^^^"*'. In d < 4 dimensions also the 
part of the pinning potential related to ps becomes rele- 
vant which gives the dominating contribution to the FL 
displacements. Both, a Flory argument Q and more de- 
tailed calculations for a related one dimensional prob- 
' 23I 24, [2I] give in d = 3 dimensions SG{x,y,z) ~ 
x\/Ld _ the related study [l^l Villain and Fernandez 
found from a non-perturbative RG that for d < 4 the 
disorder renormalizes to strong coupling. We will show 
below that this case gives qualitatively the same results. 

To get more information about a real 3-dimensional 
system we consider next the stability of this glassy phase 
with respect to point disorder by using an Imry-Ma argu- 
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26|. The energy gain from the point disorder in a 
vesion L'^-^LyLz is of the order -((i?p(w))LyL^L''-2)i/2 



L 



27| which has to be compared with the elastic energy loss 
'^{ciie'^LyLz/L'^ + T,zLy + T,yLz) ^ L"^. If one ignores 
the fluctuations of u and replaces {Rp{u)) by a constant 
one finds that point disorder is irrelevant above d — 2 
dimensions. This critical dimension is further decreased 
to zero if the fluctuations of u are taken into account by 
using Rp{u) ^ ]J-'^-'°)l'^ _ A similar argument shows the 
irrelevance of columnar disorder. This argument applies 
for Lo where Sj^/^ has developed. 

Flux creep. Next we consider the flux creep under 
the influence of a transport current parallel to the defect 
planes which creates a driving force density f = J A B/c 
perpendicular to them. J is the current density. The 
motion of the FL bundles under the influence of f occurs 
then by nucleation of critical droplets in which FLs are 
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Thus the non-linear resistivity is reduced considerably 
with respect to the case of point impurities. A similar 
consideration for the Bose glass gives n = 1 which is, as 
far as we are aware, also a new result 

To summarize the results obtained so far we remark 
that the new phase described here is characterized by 
(i) diverging elastic constants C44 and cge but a finite 
compressibility cn, (ii) a transverse Meissner effect as 
well as a resistance against shear deformation, (iii) large 
sample to sample fluctuations of the susceptibility, (iv) 
an exponential decay of the structural correlations (in 
d = 3) and (v) a creep exponent /i = 3/2. Since the 
totality of these properties is different from the Bragg 
glass or the Bose glass, we will call this new phase a 
planar glass. This phase is also different from that found 
for equally spaced defects which is incompressible [28|. 

Strong disorder. If the disorder is strong, i.e. if Lo ^ 
£d, (we ignore for the moment the point disorder) each 
defect will be completely overlapped by a FLL plane to 
gain its full energy. Integrating out the displacement field 
between two adjacent defect planes we get in d = 3 

Cll (Ul+l - UiY 



n 



LzLy 



PqVd 



^ ^ ^iGon{xi-Ui} I 



where Gd — 27r/^, Axi+i = Xi+i — Xi and the sum 
over i is over the defect planes. For w^i — > 00 we have 
Xi — Ui = iui with Hi integer to minimizes the pinning 
potential. Minimizing subsequently the elastic energ y al - 
lows the exact determination of the ground state [25|: 

— Xi — •^X]j=i [-^^j/^Ig' where [x]g denotes the clos- 
est integer to x. For £d ^ £, <5'g(i") is again decaying 
exponentially in the a;— direction on scale £d. Consider- 
ing flux creep due to a driving force f perpendicular to 
the defect planes in d = 3 we obtain then the same form 
of the non-linear resistivity Eq. ([5]) as in the case of weak 
disorder. This formula applies for small currents where 
droplets cover many planar defects. Thus both weak and 
strong disorder give the same results for the correlations 
and the flux creep. 

Displacement parallel to the defects, dislocations. Next 
we include displacements Uy parallel to the defects. In 
the case of strong disorder each defect is occupied by 
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a single FL layer and hence Ux{xi,y, z;ni) = Xi — lui, 
Vy, z, to maximize the pinning energy gain. Even without 
point disorder we obtain then a non-zero displacement 
Uy. This can be seen most easily in the isotropic case 
where ada:Ux = -dyUy, here a = (cn - C66)/(cii + cge) 
is the Poisson number, < ct < 1 [2^. The strain dxUx 
in the segment between the defects at Xi+i and Xi is 
dxUx « 1 - Mni+i/Axi+i where An^+i = (n^+i - Ui). 
The difference of the strain dyUy in neighboring segments 
is then AdyUy « a£[Ani^i/ Axi+i — Aui/Axi] which is 
of the order ztal/in- On the scale Ly this implies Auy ~ 
±a£Ly/£D- To avoid a diverging shear energy one has to 
allow for dislocations with Burgers vector parallel to the 
y-direction sitting at the defects. Their distance in the 
y-direction is of the order -^d/ct. Comparing the energy 
of an edge dislocation piercing the crystal to the energy 
gain from the disorder we find that dislocations will be 
present if ac^ei^^c <C iovo- In general, the network of 
additional FLL sheets spanned by the dislocations will 
be complicated. The resulting state is ordered in the 
sense that Sy,!]^ are non-zero and hence the transverse 
Meissner effect is still present. 

Adding weak point impurities will further randomly 
shift the positions of the dislocation leading most likely to 
decay of translational correlations in the j/z— plane. Since 
the Burgers vector of the dislocations is always parallel 
to the defects, creep in the x-direction is not facilitated. 

and Sz are still both non-zero and hence we recover 
the creep law Eq. ([5]). To describe creep parallel to the 
defects one has to take into account the interaction be- 
tween the dislocation, a situation not considered so far 
[sot . We leave this case for further studies. In the case 
of weak pinning qualitatively the same behavior can be 
expected on scales 3> Lo, in particular if the flow is 
again to the strong coupling fixed point. If the samples 
exhibits orthogonal families of (non-intersecting) defects, 
long range order in the xy— plane is destroyed even with- 
out point disorder on scales larger than Lo- The creep is 
now limited by the slowest mechanism and hence Eq. ([5]) 
is likely to be valid for all current directions in the xy- 
plane. 
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